Single-site Anderson model. II Perturbation theory of symmetric model 
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The strong electron correlations caused by Coulomb interaction of impurity electrons are taken 
into account. The infinite series of diagrams containing irreducible Green's functions are summed. 
For symmetric Anderson model we establish the antisymmetry property of the impurity Green's 
function, formulate the exact Dyson type equation for it, find the approximate correlation function 
Z a (iui) and solve the integral equation which determines the full propagator of the impurity electrons. 
Analytical continuation of the obtained Matsubara Green's function determines the retarded one 
and gives the possibility to find the spectral function of impurity electrons. The existence of two 
resonances of this function has been proved. The smooth behaviour was found near the Fermi 
surface. The two resonances situated symmetrical to the Fermi surface correspond to the energies 
of quantum transitions of the impurity electrons. The widths and heights of these resonances are 
established. 

PACS numbers: 78.30.Am, 74.72Dn, 75.30.Gw, 75.50.Ee 
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I. INTRODUCTION 

We shall investigate the properties of the normal state 
of single-site Anderson model. For that we shall use the 
results of diagrammatic theory for this model developed 
in our previous paperW. In that investigation the no- 
tion of irreducible Green's function has been introduced. 
These functions contain the main spin, charge and pair- 
ing fluctuations caused by the strong Coulomb repulsion 
of impurity electrons. We have determined the notion of 
correlation function Z aa i as composed from strong con- 
nected diagrams containing the irreducible Green's func- 
tions. In superconducting state of the system there are 
additional correlation functions Y aa i and Ys„i which are 
the order parameters of the system. The correlation func- 
tions are the main elements of the Dyson type equations 



for one-particle renormalized Green's function of con- 
ducting and impurity electrons. We shall restrict our- 
selves only by the discussion of the properties of the nor- 
mal phase of the system and determine the corresponding 
Green's functions of conduction and impurity electrons: 



(1) 



G^(k,r|k',r') = -(TC kff (r)C kV (r')i7(/?)); ; 

<w(r|r') = -(tu(t)JAt')u(P)) c , 

where index "c" means connected. 

Fourier representation is denoted as G aa i (k, k'|iw) and 
g aa i(iu)) correspondingly. For these two functions we 
have obtained the results M : 
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G CT (k,k'M 



<WG°(kM + ^G°(kM 3(T MG°(k'M, 

A g (w) 
g° a {iuj) + Z a (ioj), 



(2) 



where zero order propagators of the conduction and im- n _ _ exp( /?£/) + cxp [ @(2ef + U] 

purity electrons have the form (a = —cr). Zq 

Zq = l + 2exp(-/3e / ) + exp[-/3(2e / + [/], 

GSM*,) = (—(k))-, G o M = i y> k | 2 G°(k M = [YM£2m. 

q , . . 1 n a n a IS J tu — e 

IL) — €j IUJ — Ef — U 
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Here e(k) is the energy of conduction band and e/ of 
local impurity electrons, po(e) is the density of states 
of the bare conduction band and matrix element of hy- 
bridization Vk is supposed dependent of the energy. U 
is Coulomb repulsion of the impurity electrons, u = 
uu n = (2n+l)n/(3 is odd Matsubara frequency. The equa- 
tions (2) are exact, but for correlation function Z a (iuj) 
doesn't exist exact Dyson type equation and only the ap- 



where 

G °^) = Jf E V^G^Miu). (8) 

kik 2 

On the base of equations (2) and (3) the last function (8) 
can be presented in the form: 

(9) 



In the second Section of this paper we shall discuss 
the simplest case of symmetric impurity Anderson model 
with the condition 2ej + U = 0. In Section III the spectral 
function of the impurity electrons is analyzed and the last 
Section IV contains the conclusions. 



proximate contribution can be available: see Fig. 9 of 
paperW . Our main approximation formulated in paperW 
comes to the summation of the ladder diagrams which 
will be enough to obtain the main contributions of the 
spin and charge fluctuations. This approximation has 
used only the simplest irreducible Green's function g^ lr 
which is iterated many times. It has the form: 



(4) 



(5) 



(6) 



(7) 



By using the definition (2) of correlation function of nor- 
mal state A a (iuj) and approximation (7) for function 
Z a (iui), we obtain the final integral equation for A CT : 



(10) 

I 

II. SYMMETRIC MODEL 

In symmetric case when ej = — U/2 < and ef + U = 
U/2 > we have more simple equations: 

0/- ^, iuj 





Z aa ,{T-r) = - ^2 J dr i J ^ rfT 232 0)^r [ CT J T ; CT l> r ll CT 2,r2;c^^T']^y k * l y k2 G (J2(Tl (k2,r2|kl,Tl), 

klk2<Tl<T2 o 

or in Fourier representation 

Z aa ,{iu) = -4X]E E 92 0)^r [ CT ^ w ; CT l^ w ll CT 2,^Wl;CT^^M^K' 1 ^k2G CT 2 ( T 1 (kl,k2|^Wl). 

" ui cricr 2 kik 2 

Here we take into account the conservation law of the frequencies: 

g^ lr [a, iuj; a\, iu>i\ai,iu>i; a' , iuj'] = f35(oj — uj')g^ lr [a, iuj; a\ 1 iuji\a2 1 iuji; &' , iu]. 
In paramagnetic phase we have more simple equation (a' = a): 

Z a {iuS) = - 4 ^ ^ 92 )%r K ioJ ; ai , iu>i \ai , iui ; a, iuj]G ai ) , 



l-A CTl (i Wl )G0(i Wl ) 
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Z = 2(l + exp(-/? e/ )) = 2(l + exp(/?£//2)), 
n CT = 1/2 (11) 

and the antisymmetry property of zero order impurity 
Green's function g®(— iuj) = —g®(iui) takes place. 

Additionally we suppose also the evenness of the ma- 
trix element V(e) and of the bare density of state Po(e). 
In this case the function G a a (iuj) is also antisymmetric 
Gl{-iw) = -G»(i W ). 

Thanks these antisymmetry properties we shall look for 
the antisymmetric solution 



A a (-icj) = -A a (iu>) 



(12) 



of the equation (10). Analytical continuation of these 
functions have the property of oddness of their real parts 
and evenness of imaginary parts in conformity with equa- 
tions 

g° a (E + iS) = - g °(-E-iS), 



G° a (E + i5) = -G° a (-E-i5), 
Al(E + iS) = -Al(-E-iS), 



(13) 



where 



G° a (E + iS) 
1(E) 
T(E) 



1(E) - iT(E), 

\V(e)\ 2 Pa (e)de 
E-e ' 
7t\V(E)\ 2 Pq (E). 



(14) 



1(E) is the principal part of the integral. This func- 
tion is antisymmetric. T(E) is the band width of the 
virtual level and an even function of energy. The sym- 
metric impurity Anderson model has the advantage to be 
of a simple form for the irreducible two particles Green's 
functions in different spin and frequency channels. In 
this special case we have' 2 ' 3 ! 



g 2 [<J, iuj; a, \a, tu>i;a, iuj\ 



f3(U/2) 2 (l-5ojLJ 1 ) 



[(iu) 2 - (U/2) 2 ][(iuji) 2 - (U/2) 2 



(15) 



u I 

pU5(oj - wi) exp(/?[//2) 
" (1 + exp((3U/2))[u 2 + (U/2) 2 ][iv 2 + (U/2) 2 

+ 4(C//2) 2 



(1 - eMPU/2)) 



[a^a^a^a^ = - {2 (1 + cxp(/3c//2))[w2 + (f//2)2][w 2 + ([//2)2] 

/3US(to +wi) 



+ ■ 



[lo 2 + (U/2) 2 }[cj 2 + (U/2) 2 _ 
1 



(l + exp(/3t72))[u, 2 + (U/2) 2 ]^ 2 + (U/2) 2 ] 
1 



^ 2 + (U/2) 2 ][u 2 + (U/2) 



212 



+ 



(16) 



^ 2 + (l72) 2 ] 2 [^ 2 + ([7/2)2] 



Now we come back to equation (10) and note that thanks the antisymmetry property (12) of A a (iu) function only 
those terms of equations (15) and (16) which contain Kronecker 8 - symbols give the non zero contribution in the 
right hand part of it. The result of summation has the form: 



A a (ius) 



+ 



Z(U/2f 



(ioj) 2 -(U/2) 2 [(ioj) 2 -(U/2) 2 ] 2 [l-K a (ioj)G°(iio)Y 



(17) 



We notice that the scattering channel with opposite spins (16) gives in equation (17) the twice contribution in 
comparison with parallel spin channel (15) and both of them arc added together giving the factor 3 in the right - 
hand part of equation (17). There are two solutions of equation (17) and we take that of them which has the correct 
asymptotic behavior X,y(iuj) — > J- when \w\ tends to infinity. This solution has the form: 



A a (iu>) = 



2G0M[M 2 -(L72) 2 J 
iLjG° a (uj)]y/l-12Q(iLj)} . 



{[M 2 - (U/2) 2 + ^G°M] - [M 2 - (U/2) 2 



(18) 



r 



where 



Q(iui) 



(172)G°M 



(iu:) 2 - (U/2) 2 - iwG°(iw) 



We have used that branch of square root which gives 
one when Q(iu) tends to zero. On the base of equations 



(19) 
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(2) and (8) we obtain the renormalized impurity electron propagator 



g(iu>) = 



2[.g°M[G°M]- 1 +3QM(g°M - [G° (im)}- 1 ) 2 } 

[[GO(iw)]- 1 - s°M][l + v/T - T2QM + 6QH( 5 »HG»M-1)]' 

I 



(20) 



In the last equations the spin index a is omitted because 
it is not significant. Equation (20) has been obtained 
by taking into account the spin and charge fluctuations 
contained in the correlation function Z a (iuj). 

The spectral function of the impurity electrons is equal 

to 



Af (E) = -2Img(E + iS) , 



(21) 



where g(E + iS) with S = +0 is the analytical contin- 
uation of the Matsubara to retarded Green's function. 
In absence of the correlation function Z a (iui) instead of 
equation (20) a more simple form appears 



l-c/0 



l 



[^mp 1 - 



(22) 



which can be named as Hubbard I approximation. This 
equation contains the zero order Green's function g (itj) 
determined by equation (11) and averaged by hybridiza- 
tion conduction electron function G®(iu) of equation 
(14). 

Analytical continuation of equation (22) gives 



gi(E + iS) 



E - (U/2) 2 E~ 1 - 1(E) + iT(E) ' 



(23) 



We shall compare our renormalized spectral function (21) 
with more simple bare quantity A® (E) : 



A° f (E) = -2Img°(E + iS) = 2k5(E - {U/2?E- X ) = tt (s{E - ^) + S(E + 
and the value Aj(E) obtained in Hubbard I approximation: 



Af(E) = -21-mg 1 (E + iS) 



2T(E) 



(E - (U/2) 2 E- 1 - 1(E)) 2 + T 2 (E) 
I 



(24) 



(25) 



Here the odd 1(E) and even T(E) functions are deter- 
mined by equation (14). Two resonances of equation (24) 
situated at energies E = ±-j have not the width. After 
some interactions taken into account by Hubbard I ap- 
proximation a new spectral function (25) appears. It has 
two resonances with shifted values of energies E = ±Eq 
determined by the presence of function 1(E): 



En 



En 



I(E ) = 0. 



(26) 



These resonances are broadened by the presence of the 
function T(E), which is the width of the virtual level. 
This function T(E) determines the height and width of 
the both resonances. Near the new values of resonance 
energies ±En we can approximate (25) with more simple 



Lorentzian forms: 



2r 



where 



V> 2 (£T£o) 2 + r 2 ' 



^ = l + (^r?-l'{E Q ) 

Atjn 



(27) 



(28) 



and Eq is determined by equation (26). 

The both approximations (24) and (25) give the zero 
value of spectral functions on the Fermi surface, where 
E = 0. The correctness of the result has to be verified 
by the sum rule 



oo 

/ 



A f (E)dE = 2tt. 



(29) 
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Equation (24) fulfils this condition. If we omit for the 
simplicity function 1(E) in equation (25) we can verify 
also the fulfillment of this condition for equation (25), but 
not for its approximations (27), because the parameter 
ip is not equal to two. Some details connected with the 
choice of the density states can be found in Appendix. 

In the next Section we use the equations (19) and (20) 
to obtain more complete spectral function of impurity 
electrons and to verify the existence of the resonance at 



zero energy 



III. RENORMALIZED SPECTRAL FUNCTION 

Analytical continuation of equations (19) and (20) 
which are necessary to calculate the spectral function 
have the form 



g(E + iS) = [[(^(E + iS^-g^E + iS)]- 1 x 

2 [g°(E + i8)[G°(E + iS)}- 1 + 3Q(E + iS)(g°(E + iS) - [G°(E + iS)}- 1 ) 2 ] 



1 + sf\ - 12Q(E + iS) + 6Q(E + iS)(g°(E + iS)G°(E + iS) - 1) 



(30) 



UG°(E + i5)/2 \ and supposing the smallncss of the parameter 2T/U we 



Q(E + iS) y^ 2 — (U/2) 2 — EG°(E + iS) ) ' can approximate equation (30) by more simple one: 

First of all we shall analyze the behavior of the renor- 
malized Green's function on the Fermi surface for E = 0. 
Near the value of the energy E = we can approximate 
the quantity Q(E + iS) by the expression 

Q(E + i5)\ E = =-(?^) (32) 



g(E + i6) ~ 



For little values of energy we have 



6Q(E + iS) 



G°{E + iS)[l + y/1 - 12Q(E + iS) - 6Q(E + i6)} ' 



g(E + iS) ~ 



6(2r/[/) 2 



T 2 {E) +E 2 I'(0) 2 



1 + a/1 + 12 (2T/U) 2 + 6 {2T/U) 2 



-EI'(0)+i(T(E) + -^-) 



(33) 



(34) 



The spectral function of impurity electrons in the region of little values of energy E is different from zero and has 
the form 



ME) 



i2(2r/io 2 /r(o) 



1 + 6 (2T/U) 2 + y / i + i2(2r/c/) i 



(35) 



In Appendix are cited the values of the quantity I'(0) 
equal to 4p(0)V 2 (0)/W or nV 2 (0) p(0) / D for two differ- 
ent model of density of states. The result (35) differs 
essentially from the zero value obtained in such more 
simple approximations for Af(E) as A®(E) and Aj(E). 
Thus we have established the existence of two resonances 



at energies E = ±Eq determined by (26) and the pecu- 
liarity at E = 0. We can find the corrections for spectral 
function Af(E) for two resonances E = ±E - With this 
end in view we determine the quantity Q for E = E a 

Q[E + iS) = - (I(E ) - iT(E )f ( 2Eo r {Eo) ) , (36) 
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where we suppose that the value Eq is inside the edges 
of the band width. After some calculations made 
with the supposition that quantities I(Eq) /T(Eq) and 
UI{Eq)/2EqT{Eq) are large and are superior in number 
to one we obtain 



1 



4£n 



V3(1 + 3(U/2E ) 2 ) UI(EoY 



(37) 



This quantity essentially differs from the value Ai(Eo) = 
2/T(Eq) being considerable less than it. Thus the process 
of renormalization of the impurity electron propagator 
results in appearance of the peculiarity at E = and to 
diminution of two resonances situated at E = ±Eq. 

The full spectral function can can be calculated on the 
base of the equation (30) and (31). On the Figure Q] 
the results of numerical investigation of the full spectral 
function are presented for different values of theory pa- 
rameters. For comparison on this figure also the result 
obtained in Hubbard I approximation is presented. As 
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can be seen from Figure [T] there are two sharp resonance 
peaks near the energies E = ±E and smooth behaviour 
near E=0. The distance between peaks is determined 
by parameter U and the height and width of peaks are 
determine by parameter T. 

IV. CONCLUSIONS 



We discuss the symmetric Anderson impurity model 
and take into account the strong electronic correlations 
of the impurity electrons by elaborating the suitable dia- 
gram technique. This paper is based on the previous our 
oneW founded on the diagrammatical investigation and 
analysis of the properties of single-site Anderson model. 
First of all we establish the antisymmetry property of the 
functions g°(iuj), A a (iu) and g a (iu), and use the exact 
Dyson type equations for our functions. The special ap- 
proximation for the correlation function Z a (ioj) has been 
obtained which gives the possibility to close the system of 
equations and to find the solution for renormalized func- 
tion g a (iu!). This Matsubara Green's function is contin- 
ued analytically to obtain the retarded one. 

Spectral function of impurity electrons is obtained and 
the structure of resonances and their properties are an- 
alyzed. Two of resonances of this function at E = ±Eo 
(Figure [1]) correspond to the energies of quantum transi- 
tions of single - site impurity and the smooth behaviour 
was found at the energy E — 0. The details of the spec- 
tral function renormalization are based on the proper- 
ties of real 1(E) and imaginary T(E) parts of function 
G^(iw), which is conduction band electron function av- 
eraged by the hybridization interaction. The values of 
these functions and of the values of energies ±Eo are 
discussed in a special Appendix. 
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FIG. 1: Spectral function Af(E) x W for different values 
of the theory parameters as function of energy E/W in the 
Hubbard I approximation (case (a)) and in the our ladder 
approximation (cases (b) and (c)). In the cases (a) and (b) 
U = leV, W = 2eV, V = 0.1 eV and in case case (c) 
U = 1.5eV, W = 2eV, r = 0.15eV. 
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APPENDIX A: SIMPLE EXAMPLES OF 
DENSITY OF STATES 

We can demonstrate some simple examples of the 
choice of the density of states and of the correspond- 
ing functions 1(E) and T(E). For simplicity the energy 
dependence of the matrix element of hybridization V(e) 
is supposed smooth and is neglected. One example of 
density of states has been proposed in the paper' 4 !. In 
this paper the following equations are used 

po(e) = Po(0) (l - (e/uf) , |e| < W, 

1(e) = p (0)V 2 (0)[2e/W+(e 2 /W 2 - l) In 

r(c) - nV 2 (e)p (e). 



W 

~w 

(Al) 



where 2W is the conduction band width. For little value 
of energy we have /(e) = l'(0)e with 



J'(0) = 4 Po (0)V 2 (0)/W 



(A2) 



and for E — > ±00 function 1(E) tends to zero as 
1/E. In the case (Al) the equation (26) takes the form 
(x = E /W): 



x - a 2 /x a - bip(x ) = 0, 



(A3) 



where 



(p(x) 



2x + (x 2 - 1) In 
p(0)V 2 /W. 



x + 1 



; a = U/2W; 



We note that the functions po(e) and T(e) exist only in- 
side the edges of the conduction electron band \E\ < W 



whereas the function 1(e) and the solution of equation 
(A2) can exist also for \E\ > W. Therefore we have to 
find the solution of (A2) with x < 1 and consider the con- 
ditions for the values of parameters a and b compatible 
with this requirement. 

Another simple example of the density of state is one 
with Lorentzian shape' 5 ' 



p (e) =2D/(E 2 + D 2 ) 1 



(A4) 



with chemical potential placed at the e = 0. This choice 
has the advantage of not introducing band edges. It has 
the parameter D as an effective band width: 



1(E) = 2irV 2 (0)E/(E 2 + D 2 ) 



T(E) = tt P o(0)V 2 (0) = 2irV 2 (0)/D. 



(A5) 



In this case instead of equation (A3) we have other values 
of parameters a and b and other form of function ip(x): 



x-a 2 /x-bx/(l + x 2 ) = 0, 

a = U/2D; b = 2nV 2 (0)/D 2 . 
Equation (AQ) has two solutions ±xq with 



(A6) 



X 



V(l - a 2 + b) 2 + Aa 2 + a 2 + b - 1 /y/2 



1/2 



(A7) 

The value of parameter I'(0) is equal to 2ttV 2 (0)/D 2 = 
7tV 2 (0) Po (0)/D. 



